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ABSTRACT: We consider the spin determination of new colored particles in the miss- 
ing energy plus jets channel at the early stage of the discovery. We use a three 
site moose model to describe the low energy Lagrangian of all same spin partner 
(LHT or UED like) models and check the gauge invariance of the amplitude. For 
the benchmark production and decay channel pp — U“U® — uuBy Bu, in con- 
trast to those in supersymmetric models, there are spin correlations which affect the 
polar and azimuthal angle distributions of the quarks from the heavy partner U® 
decay. We show such effects would be visible in the E-pmiss/Meg distribution and the 
reconstructed azimuthal angle correlation using MAOS reconstruction. 
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1. Introduction 


The existence of dark matter in our Universe is solid evidence of new physics beyond 
the standard model. Dark matter must consist of a new elementary particle, whose 
lifetime is much longer than the age of our Universe. One way to incorporate the 
dark matter candidate beyond the Standard Model (SM) is to introduce a parity 
structure. This parity structure can be incorporated in well-studied models such as 
supersymmetry (SUSY) (R-parity), little Higgs (LHT) (T-parity) [1], and universal 
extra dimensions (UED) (Kaluza Klein parity) B, B, Ø. Under such a parity trans- 
formation, all the SM particles have even parity, their partners have odd parity, and 
the interaction vertices multiplicatively conserve the parity. The lightest parity odd 
particle becomes stable, and becomes invisible in collider searches. 

Observation of large missing transverse momentum at a collider would raise fur- 
ther questions about the properties of the underlying physics. The signatures of 
SUSY and same spin-partner models are quite similar ff, 6]. In supersymmetric 
models, scalar quarks and scalar leptons have spin 0 while in LHT and UED, quark 
and lepton partners have spin 1/2. The gauge boson partners are spin 1/2 in super- 
symmetric models, while they are heavy gauge bosons in LHT and UED. Therefore, 
it is critical to determine the spin of partners of the SM fields in order to distinguish 
different types of theories. To distinguish the spin structure of SM partners, we need 
to study a quantity that can be measured in the early stage of the Large Hadron 
Collider (LHC). At a hadron collider, the main production channels contain colored 
objects, whose main decay channels involve jets. The total cross section may serve as 
an initial hint to the spin of the odd partners [/]. However, at this stage of discovery, 


small statistics make it impossible to reconstruct the mass, which the cross section 
depends on strongly. 

The spin-dependence of decay patterns of new particles is useful in determining 
the spin structure. The decay distributions depend on the spins of the initial, interme- 
diate and final state particles. For example, in the two leptons + jets + ETmiss chan- 
nel at the LHC, the spin dependence appears in the jet invariant mass distributions. 
The simple my distribution of dP /my o my arising from a Gg > x$j > iil > x ofl 
cascade decay in SUSY reflects the fact that the slepton is a scalar. In addition, the 
charge asymmetry in the m; distribution arises from the polarization [9] of the 
spin 1/2 neutralinos. The polarization arises because the LHC is a proton-proton 
collider, so squarks are produced more often than anti squarks. 

The difference between the invariant mass distributions of SUSY signatures and 
those of same spin partner models have been studied previously 14]. 
The spin and interaction dependence is sizable, but such studies at the LHC may not 
be straightforward, because of the difficulty of selecting the correct jets and leptons. 
For a relatively clean jll signature, the decay branching ratio of g — x$ is typically 
only around 30%, and y2 — 1 — x? is typically less than 10%. It is not the dominant 
decay channel and is therefore statistically limited. To distinguish the spin structure 
of quark partners in the relatively early running of LHC, we need to focus on the 
polarization dependence of the dominant decay channel. 

The other useful information lies in the spin correlation of particles from different 
decay chains. Here invariant mass distributions are not the useful quantity, rather 
the distributions of the polar angle and the azimuthal angle correlations [15 (16, 
of decays. Due to the uncertainty of initial parton energy, it is not easy to see 
the angular correlation exactly at the LHC, but the effect should exist in various 
distributions in cases with strong angular correlations. 

In this paper, we study the spin correlations of the qq — QQ process where 
Q is a spin 1/2 quark partner. We use a three-site deconstructed model which 
universally parametrizes low-energy effects in popular same spin partner (LHT or 
UED like) models with gauge boson and fermion partners. With such a simplified 
gauge symmetry breaking sector, it is straightforward to check the gauge invariance 
of the amplitude which guarantees good high energy behavior. In our model, the pair 
production and the decay of the spin 1/2 quark partner can have a spin correlation 
which is in contrast to the isotropic decays of scalar quarks. Among the possible 
signatures, the two high pr jet + missing Er +X channel is promising, where X 
is an arbitrary number of jets and leptons. The quark partners can directly decay 
into a light stable gauge boson partner and a quark, and the two high pr jets come 
from their direct decays. It is essential that the jets coming from the two body decay 
tend to have a large pr because it makes their selection much easier. We stress that 
the channel we are proposing, which is one of the main discovery channels at the 
LHC, could be the most promising one to see spin correlations. Once the integrated 


luminosity becomes closer to 1 fb~!, the LHC will be able to access partners with a 
production cross section above O(1) pb~+. Therefore we propose simple distributions 
that are sensitive to the spin correlations, so that we would be able to determine the 
spin structure as early as possible. 

The paper is organized as follows. In Section 2], we introduce a three site model 
which gives us the low energy spectra of the massive gauge boson and fermion part- 
ners without model parameter constraints. We double check the gauge invariance of 
the amplitudes which guarantees their good high energy behavior. In Section f, we 
use two sample spectra in parton level simulations and present the polar angle and 
azimuthal angle dependences. We also calculate the amplitudes in different helicity 
channels. In Section ff, we demonstrate how the jet level distributions are affected 
by the spin correlations. Section [| contains our conclusions. 


2. Quark partner production in a three-site moose model. 


There are various models which contain SM “partners” with the same spin and an 
odd Zə parity. Examples such as the universal extra dimension model ( “bosonic 
supersymmetry” ) [P| B, ff], its variants with split mass spectra [9], little Higgs 
models with T-parities ||| or even the warped models with KK parities belong to 
this group 20]. In the process pp > uru: followed by u; decay into a quark and 
light stable particle (+ softer jets or leptons), one can safely neglect the effects 
of electroweak symmetry breaking, and all the massive particles gain their masses 
through the extra gauge symmetry breaking or the compactification with bulk gauge 
symmetry. From the dimensional deconstruction [B], point of view, all of these 
models can be described by the three site moose diagram. Another reason we choose 
the three site moose digram here is because we need to describe a Goldstone boson 
with relatively simple mass mixings, which are needed in the Re gauge to check the 
gauge invariance of the amplitude and the high energy behavior. On the other hand, 
it would be very cumbersome if we were to use the the fifth component of the gauge 
boson in the language of extra dimensions. 

We consider a three site U(1) moose model which mimics the low energy theory 
for the scattering wu — uu; in any kind of same spin partner (LHT or UED like) 
model. For the non-abelian case, our result would only differ by a color factor. The 
moose diagram is presented in Fig [I]. The gauge couplings in each moose are gp’, ga 
and gp, and we set gp = gp so the model has a Zə parity. The fermions charged under 
the gauge groups U(1)g and U(1)p have equal Dirac mass terms —Mw?w2 + h.c. 
and -MYE pe + h.c.. The two bifundamental scalar fields © and X’, which we 
will call the “link” fields, are charged under gauge groups U(1), U(1)4 and U(1)a, 
U(1)p with charge (1,—1) and (—1,1) respectively. They couple to the fermions 
through the Yukawa couplings with the same coupling strength yarns + h.c. and 
ype dye + h.c.. 
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Figure 1: The moose diagrams for the three site split-UED model with a single 5D Dirac 
Fermion. An arrow into a site means that the particle transforms under the fundamental 
representation of the relevant site and an arrow out of a site means that the particle 
transforms under the anti-fundamental representation. The solid lines stand for Weyl 
fermions and the dashed lines represent scalars. The left panel shows how to get a left- 
handed chiral zero mode and the lower right shows how to get a right-handed chiral zero 
mode. 


The link fields get their vacuum expectation values (2p) = (Xiz) = udjg and 
spontaneously break U(1),, U(1)4 and U(1)z into the diagonal group U(1)9. The 
kinetic term for the link fields Tr[(D,,©)'(D,,©)] + Tr[(D,»’)'(D,,5’)] generates the 
mass terms for the massive gauge bosons. The mass matrix of the gauge bosons is 


gpu?  —gagpu’ O By 
(B, A, Bi) | -gagsu? 29%u? —gagsu? | | A, | - (2.1) 
0 —gagpu ghu? Bi 


We can find the mass eigenstates of the gauge bosons in terms of the gauge 
eigenstates 


G, sg/ V2 Cg sgi V2 B, 
G! | =| -1/72 0 1/72 | |4,]- (2.2) 
G? Cy/V2 —8q c,/V2 Bi 


The inverse transformation between the gauge eigenstates and mass eigenstates is 


Bu sg/ V2 -1/V2 cg/ V2 G, 
A, | = Cg 0 =s; GI] à (2.3) 
Bi} \sy/VB 1/vE oI NE 


where s, = V2ga/9, Cy = gp/g and g = \/29, + g2. The corresponding gauge boson 
masses for Gi Gy G? are 0, ggu and gu. 


Similarly, the gauge invariant Dirac mass term and the Yukawa interactions give 
the fermion masses. The fermion mass matrix is 


Mr yru O\ (YÈ 
(vp vi ve) { 0 0 0 ve (2.4) 
0 yru My, a 
We can also find the mass eigenstates of the fermions in terms of the gauge eigenstates 


Y? s¢/V2 —cy sp/ VN (vP 
yt | = | 1/2 0 12| (ee 1 (2.5) 
wi cs/V2 sg ce/V2) NWP 


Wr O 1 0 YR 
p=- Aa] | vd | . (2.6) 
v2) \ av 01/v2) Wë 


The inverse transformation between the gauge eigenstates and mass eigenstates is 


L 85/V2 —-1/V2 cr/V2\ [V 
L = =Cf 0 Sf pi , (2.7) 
P sp/V2 1/v2 erva Ni 


ve 0 —1/V2 1/ VÀ (ve, 
wA |= [1 0 0 bp |. (2.8) 
UR 0 1/V2 1/2) WW} 


Here, the fermion mixing angle is defined as sp = V2yu/./M? + 2(yu)? and cs =M/ 


vos 


M? + 2(yu)?, and the corresponding fermion masses for w°, Yt, w? are 0, M, 
M? + 2(yu)?. 

In this case, the model contains one gauge boson sector and one chiral fermion 
sector. The gauge boson sector contains one parity even massless state (0 mode), 
one parity odd massive state (1st mode) and one parity even massive state (2nd 
mode); The chiral fermion sector contains one parity even left-handed chiral fermion 
(0 mode), one lighter parity odd massive state (1st mode) and one lighter parity 
odd massive state (2nd mode). In order to mimic the SM and cancel the anomaly, 
one can also introduce another fermion sector which has a parity even right-handed 
chiral fermion as in the right panel of Fig. [I]. 

For the Goldstones, the interactions with fermions are 


= ypf Eyr + ybPo'YR + hc. 
C C T 
a CE Eph + hee. 


> a YR- 4/2 


cry cry Cry -CFY 
5 iva — ttn — oh! + Zhe 


v2 v2 v2 v2 


D —icpyh YRT + icpyhRYLT , (2.9) 


where the symmetry breaking of o and o’ is linear 


u=(uto)+iz. 
X = (u+) + in’ (2.10) 


and Gi, eats the linear combination 7 = (a! — m) //2. 

Now we can derive the Feyman rules for the model and calculate the full ampli- 
tude. The relevant Feynman rules that we use in the calculation are presented in Fig 
Ø. The wt — y? — G! coupling is giz = gcgs¢/V2, and the Goldstone yh — y} — m 
coupling is gip = —cry. For the scattering process wu — u ui, one must make sure 
that the calculated amplitude does not have bad high energy behavior. Indeed, as 
t Hooft and Veltman in their proof of the renormalizability of the electroweak theory 
[23], it is the gauge invariance of the theory which guarantees the cancellation of bad 
high energy behavior of the amplitude 24]. In order to check the gauge invariance of 
the amplitude, we adopt the general Re gauge and add the full amplitude including 
both massive gauge boson and Goldstone boson exchange. 

The amplitude coming from the t-channel exchange of the massive gauge boson 
in the Re gauge is 


iM = (t I) a p')Yu Pru Dz x (—) 
(o - 7 E 0-8) alk Prulk) . (2.11) 


We can extract the gauge independent piece from the gauge boson propagator as 


i v oe) a (<5) 
a Be, 2.12 
q? — m2, ( m2 P — m2€ \ mi, ( ) 


So the additional gauge dependent piece can be written as 


1 —i -m 


* == 2 fal P 1 P 
iMg 5 (9854) u(y’) HOD) ee oa mee mz, u(k’)Piu(k) 
1 
= (ycp)’U(p') Pru(p)>—zu(k') Pru(k) . (2.13) 
q= mE 
Notice that mı = M = V2yu/ts and mg, = ggu. 
The amplitude for the Goldstone interaction is 
’ _ i 
iMg = (moop ap Peu ae A=] 
u(k’)Pru(k) , (2.14) 


which cancels the contribution from the gauge dependent piece from the t-channel 
gauge boson exchange. 


Figure 2: Feynman rules for the interactions used in calculating uu > u1u1. The mixed 
coupling strength for the gauge boson and fermions is gi, = gcg8q/ V2 and the mixed 
coupling strength for the Goldstone and fermions is gj, = —ycy. 


3. Simplified Lagrangian and parton level distributions 


With the above three site model, we can adjust the input parameters MŁ, 9A, 9B, YLU, 
Yru, Mr to get the free mass spectra for the massive gauge bosons and left-handed 
or right-handed chiral fermions and their massive vector partners. The relevant 
particles in same spin partner (LHT or UED like) models with the two jets plus 
missing energy signature are 


(1) a massive gluon partner Gy or electroweak boson partner Wy. 
(2) a light stable gauge boson By. 


(3) light fermion partners QM), Q®, L™, L. (Here the superscript indicates the 
whether the fermion partner is partner of left-handed or the right-handed SM 
fermions, not the chirality of the partner.). 


All the fermion partner masses, gauge boson partner masses, and their interactions 
are free, but we assume the lightest partner is By, the partner of the SM U(1)y 
gauge boson, and all partners have parity —1, so that By is stable. The other 
massive gauge boson partners, on the other hand, could be heavier or lighter than 
the massive fermion partners, depending on which model parameters we choose. The 
interactions involving gauge boson and fermion partners are 


lige = i da [oC H OPT “gn + gRGH OPT yaR + GW OM Ty, Prat 
+91, BEY LQ Yat + RBLY RO? y qr + (Lepton part) + h.c (3.1) 


Here g3r/r, gor, and gız/r stand for the couplings between the SU(3)e, SU(2)r, 
U(1)y gauge boson partners and the SM left-handed or right-handed fermions and 
their partners. For simplicity, we take them to be the same as those of the SM gauge 
bosons in the following. We do not write the Goldstone interactions which ensure 
the gauge invariance of the amplitude. 


Since the parton distribution functions of the first generation quarks are much 
harder than those of anti-quarks or gluons for protons, the production cross sec- 
tion of quark partners is much larger than that of antiquark partners. Gluon part- 
ners can decay into quark partners or anti-quark partners, however their produc- 
tion cross section is still subdominant. Therefore, light quark partner production 
pp > QUQM/Q”Q™ is the dominant production channel at the LHC. Notice 
that the subdominant production of GyGy, GyQ and QQ may involve second 
level massive gauge boson exchanges in the s or t-channel unlike QQ production. 
The cross section for those subdominant process may be significantly different from 
that of our simplified model. 

The production and decay distributions of fermionic partners are quite different 
from those of scalar quarks in supersymmetric models due to spin correlations. Here 
we demonstrate how those spin correlations appear in different helicity amplitudes. 
The amplitude for u*(p;, hi)u® (pi, hi) — UH (pr, h U (p'r, hp) would be written 
as follows 


iM (uu —> UYR) = ig? Y2 aðan — y x 
| oe MBy 
tin; (Pp) Y Prun (Di)tin, (Pr) WPRun, (D5) 
+gluon partner exchange contributions 


+cross diagram (3.2) 


where p’s are momenta, h’s are helicities of the particles, and a,b are color indices. 
For the massive gauge boson propagator, we only write down the gauge independent 
piece since the rest of the pieces will be cancelled by the Goldstone boson exchange. 
We only show the ¢t-channel By exchange contribution and the g,, and qq” terms 
in the By propagator. The extension to the full amplitude is straightforward and 
for the numerical calculation we use the full amplitude. The amplitude arising from 
the guv term in the center of mass (CM) frame is 


E? rad 
iM, = ig? Ve barbavSn,.395.0 a ppd [26,4 Br — (=r Pp) 5 (3.3) 
— m3, 


which is non-zero for h; = h; = 1/2 and Ay = hy — h = 0. On the other hand, the 
gq’ /m term after using the equation of motion becomes 


2 
MYR) 


. - 72x 72 = = 1 1 
= —ig Y,“ba beh ——— P (pi , P (p;), (3.4 
iM2 1g Fy Cab Oa'b (2 -— ma, ym, a) R Un; (Di) ün (Py) R unt (p ) (3.4) 


which is enhanced by a factor of m7,,)/m,- Note that this coupling is Yukawa- 
like, and it therefore flips the helicity of the fermions in the relativistic limit. The 
amplitude in the CM frame is 


2 M? Eom \* 
iM, = —ig' Y dabðarv Ön, 94,09A;,0 — Z (=) 
Mp, (4 = Mpy) 


hetl 
| (25y,,-3 Br +a- BA) 
3 iO MQR) 
— 3 Dua (cos Oe (iA) (65,00 — Bs) + 2v/2|X,| r + 251 30,0 31) f 
(3.5) 
In the relativistic limit, 3 + 1, the amplitude is dominated by (hy, h) = (1/2, 1/2) 
and (—1/2, —1/2) contributions. When moe) /mg > 1, the h = —1/2 amplitude is 
dominant due to the Yukawa nature of the coupling. 


The decay of polarized Q® (Q® — qBy) is not isotropic. The amplitude 
is dominated by the h = 0 component of By if mom > Mpy, therefore the decay 


distribution is essentially that of a spin 1/2 particle decaying into a light spin 1/2 par- 
ticle and a scalar. Because of the chiral nature of the vertices such as ByUry,Prur, 
the quark in the final state is right handed. For a quark partner with h = 1/2, the 
decay amplitude is given by 


s 1+ 1 
iMp(U™ > uBy) « ee’ Tr |y,Pr pii PR Aeh + Mya) ) tAn 
2kg -pm 
_ 272 = U(R) (Ep = (3.6) 
mB 
where M is the polarization vector with n? = —1 and kj; = =n ; kg and Ep is the 


energy of By in the Q“)’s rest frame. Quarks from Q® decays with hg = 1/2 
tend to go in the direction of Q®™® helicity as a result of helicity conservation. The 
distribution is proportional to 1 + cos 6g, for hg = 1/2 in the massless limit, where 
Q4 is the polar angle of the u momentum in the U (R) rest frame, where the z axis 
coincides with the U® momentum in the CM frame. The decay distribution of QM 
(QM) — qBy or qWy) with h = —1/2 is the same as that of Q®, dr /dcosOgq x 
1+ cos @gq in the massless limit. On the other hand, it is opposite for an antiparticle 
decay. Because the number of produced quark partners is greater than the number 
of anti-quark partners at the LHC, the net polarization effect still remains in the 
signature. 

We calculate the amplitude using Madgraph/MadEvent [25] assuming that the 
collision energy for the LHC is 7 TeV. The MAOS momenta parameter cards of 
the model are obtained by BRIDGE [6]. Fully spin correlated event distributions 
are generated by calculating the amplitude at the level of the final state partons 
qq > QQ > qqBxBx. Event distributions without spin correlations are generated 
for comparison by interfacing the qq > QQ amplitude of Madgraph/Madevent to 
PYTHIA 7]. We interface the events into AcerDet for detector simulation, but 
for jet reconstruction we adopt the Cambridge/Aahen algorithm using FASTJET 
[29]. The number of events is O(10*) before the cuts. For a typical production rate 
~ pb at 1 fbt integrated luminosity, the real number of events in the data sample 
would be about ten times smaller. Nevertheless, the qualitative features that we 
show in the plots would be the same. 


In our simulation, we consider two sample points in the parameter space. We 
first choose the following spectra: The heavy gauge boson partner Gy and fermion 
partner Q™® have quasi-degenerate masses (we choose Q‘) to be relatively lighter) 
and are much heavier than the lightest gauge boson partner By. In this case, one can 
easily identify the hard jets from a heavy Q®™® decaying into By. We consider two 
different values of the light By mass, which has a large impact on the polarization of 
jets, as explained in the introduction, and illustrate how those different polarizations 
give us various different distributions. 


(A) mym = 600 GeV, mg, = 700 GeV, mp, = 100 GeV. In this case, the sum 
of the production cross section of pp > Q?QM®, QMQ and QO is 
5.8 pb. The dominant production channel Q™®Q® has a cross section of 5 pb. 


(B) Mym) = 600 GeV, mg, = 700 GeV, but the lightest gauge boson mass is 
heavier, mg, = 200 GeV. In this case, the production cross section is 2.7 pb, 
which shows that the By exchange contribution is important at Point A. 


For the spectra we choose here, the contribution of Gy exchange is larger than that 
of By exchange at point B, while the By and Gy exchange contributions are of the 
same order at the point A. 

In Fig. B|, we show the cos gq distribution for pp + U (HUA) + uuByBy. We 
generate 10+ events using Madgraph, so each histogram contains 2 x 104 entries. At 
Point A, the mass of By in the t-channel exchange is small, therefore the Yukawa- 
like coupling coming from the q“q’/ men term in the By propagator is enhanced, so 
there are equal order (1/2, 1/2) and (-1/2, -1/2) helicity states for UU), which 
tends to decrease the overall polar angle dependence. In the Point B with larger By 
mass, the Yukawa coupling decreases. Only (1/2, 1/2) helicity states for UUM 
dominate and the overall polar angle dependence is stronger. 

We now discuss the azimuthal angle correlation in the uu > U®U® = 
uuBy By distribution [6]. The amplitude is expressed as the product of the 
production and decay amplitudes, 


iM (uu > U®U® — uwuBuBu) = ip Ponai (Eom 9) Din(O*, 6") Durg] 
hyh! hh! 

(3.7) 
where P, whi = Thaw T} u and Th,w is the helicity amplitude of uu > UM Rey). 
The decay matrix D,;,(6*, ¢*) for U®) — uBy is a function of 6* and ¢*. Here 6* 
is the polar angle of the momentum of u in the rest frame of U“), where the z axis 
coincides with the momentum direction of U® in the CM frame (the production 
frame), and ¢* is the azimuthal angle relative to the production plane defined by 
the u and U“) momentum directions. The decay matrix has simple azimuthal-angle 
dependence, 


Dyn (8", $*) = Dy (O EO (3.8) 
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Figure 3: The cos 4g, distributions at Point A with mg, = 100 GeV (left) and at Point 
B with mg, = 200 GeV (right). mg, = 700 GeV and mg = 600 GeV. All histograms 
correspond to 104 events. 


The factor of exp((h — h)¢*) is 1 if the production amplitude is dominated by a 
single helicity state h = 1/2 or h = —1/2, but if both h = +1/2 amplitudes are 
sizable, the azimuthal angle correlation is important due to the large interference 
term. In that case the distribution has the form œ 1 + acos(¢) where ¢ is the 
difference of the two azimuthal angles. In Fig. H, we show the distributions of ¢ at 
Point A (mgp, = 100 GeV) and at Point B (mg, = 200 GeV). The azimuthal angle 
correlation is stronger at Point A because the amplitudes for (h, h’) = (1/2, 1/2) and 
(h,h') = (—1/2,—-1/2) are roughly equal at this point. The polar and azimuthal 
angle dependences are complementary to each other. 


4. Detecting spin correlations in jet level distributions 


In this section, we discuss the effects of spin correlations in jet level distributions. 
Unlike at ete~ colliders, at hadron colliders it is not straightforward to see the effects 
by reconstructing the events exactly, because the collision energy of each event cannot 
be measured, and the lightest gauge boson partners are missed by the detector . In 
this section, we point out a few distributions which are sensitive to the effects. 

In Fig. f], we show the distributions of Epmiss/Meg for the pp + U (RUC) process. 
In the left figure we show the distribution at point A; here the U® decay is isotropic 
(mp, = 100 GeV), which mimics SUSY squark production and decay. Fig. E] (middle) 
is the fully spin correlated distribution for pp > uuByB), at Point A, and Fig. 5] 
(right) is the distribution at Point B (mgp = 200 GeV). Here Meg is defined as 


Mer = ETmiss + X pri , (4.1) 
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Figure 4: The distributions of azimuthal angle difference ¢ of two UP) decays in the 
events at Point A (left) and Point B (right). All histograms correspond to 10+ events. 
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Figure 5: Distributions of Epmiss/Meg for pp > UUR) — urnuRBy By. Left; isotropic 
decays using PYTHIA at Point A. Middle; spin correlated decays at Point A. Right; spin 
correlated decay at Point B. All histograms correspond to 3 x 10t events before the cuts. 


where the sum is over jets with pr > 50 GeV. Meg is correlated with the typical CM 
energy of the collisions. The missing transverse energy ErTmiss is defined as 


Ermis = y (PR) + (PR), (4.2) 
where — XY) is the sum of the transverse momenta of the reconstructed objects. Up 
to the detector smearing and acceptance, it corresponds to the sum of the momenta 
of the invisible particles. 

The distribution of the events without spin correlation, (which is obtained by 
interfacing the pp — UU events to PYTHIA for isotropic decays) in Fig. $ (left) is 
clearly different from the spin correlated decay. The spin correlated events tend to 
have smaller missing energy for fixed Mag. The typical decay topologies at Points 
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Figure 6: Schematical description of event topologies of the process U®U®) — 
urnuRByBy. Left: A case of polarized quark partners (Point B). The u quarks tend 


(BVs go. Right: the case with strong azimuthal 


go in the same directions as their parent U 
angle correlation(Point A). Two u quarks tend to go in back to back directions in the 


production frame. 


A and B are shown in Fig. E] (right) and Fig. (left), respectively. The left figure 
shows a decay topology at Point B where the average polarization of the U® is 
right handed. The quarks tend to go in the directions of the parent U“) momenta, 
therefore the energy of u (By) from U“) decay is boosted up (down) in average. 
The average ETmiss/ Meg, which is correlated with P$, /my ix, is smaller compared 
with that of isotropic decay. 

In Fig. f (right), we show the other decay pattern reducing the average ErTmiss/ Mef, 
which is represented by the distribution at Point A. At this point, the average polar- 
ization of U“) is small, but the azimuthal angle correlation is stronger. In addition, 
the direction of U“) production is peaked forwards because Mp,,/Mye) is small. 
Therefore the azimuthal angle distribution of u in the lab frame and that in the 
production frame are similar. The transverse momenta of the u quarks tend to be 
back to back in the production frame, therefore the two By’s tend to be back to back 
both in the production plane and in the lab transverse plane. The By transverse 
momenta cancel on average, so the total missing pr becomes smaller. 

Next we discuss reconstruction of the azimuthal angle correlations. The simplest 
way to see this for our process is to look into the @ distribution of the leading two jets 
in the Lab frame. However, in this paper we use MAOS momentum reconstruction 
B9, to reconstruct the @ distribution. This is a method for estimating the direc- 
tions of missing particle momenta using the Mrz parameter [B2], which we explain 
below. 

The definition of Mrz is expressed as follows B2, 

Mr= min — [max (mr(pli), pT, Mes) MrR PEM) |- (43) 


pr +p5 =PTmiss 


(i 


vis 


Here pi; are the momenta of the two visible objects in the events. We adopt the 


two highest pr jets in each event, and pı, are those of the two invisible objects as 
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pe, and denote the calculated Mrz by mrz or mr2(2 jet). In this paper, we adopt 
minimal cuts to the jet rapidity 7 and transverse momentum pr so that |n;| < 2.5, 
pr, > 100GeV, and pr > 50GeV. 

The minimization is taken over all pı and py that satisfy p? + p? = pZ, and 


(p1)? = (p2)? = We 


We can define p\ and p from the p?(min) and p? (min) which minimize the 


above equation. The conditions to solve pe, are (In our case, Qı and Qo are the 
intermediate massive heavy quark partners) 


2 
(pi)? = (p2)? = Meost (4.4) 
for the transverse momentum 


pO” = pi (min), pR” = pz (min). (4.5) 
For the events with mrz ~ m™3*, pT (min) and p? (min) are very close to the true 
By transverse momenta. Therefore we define the reconstructed azimuthal angle 
difference ¢(rec) based on Po, and Pg,, namely Prec is the azimuthal angle difference 
in the CM frame of Pg, + Po, where the z-axis is along the direction of the momenta 
Pg. There are two undetermined components in the pi, and we have the two mass- 
shell equations for both sides of the decays, which give 2 x 2 = 4 solutions for 
(P, p2). Two of them give trivial solutions with (rec) = 0 , and we study the 
other two non-trivial solutions, which lead to two azimuthal angle solutions, ¢(rec, 1) 
and (rec, 2). 

The merit of using MAOS reconstruction is limited. In Fig. [f we plot the differ- 
ence between the reconstructed (rec, i) of the parton level event and the true parton 
level azimuthal angle difference ¢(true) for the events near the end point with mr > 
500 GeV. The plot is based on 20000 generated events where ~ 2500 passed through 
the cut. Here the solid histogram is A@¢(best) = min;— 2 |(ġ(rec, i) — @(true))|, the 
dark dotted histogram contains both of A¢(i) = |@(rec, 7) — d(true)| (scaled by a fac- 
tor of 0.5), and the thin dotted histogram is A¢(Lab) = ||¢,—@2|—(true)|, where ¢; 
is the azimuthal angle in the Lab frame. The matching between ¢(true) and (best) 
is good. However, once we take into account both of the solutions, ¢(rec,1) and 
@(rec, 2), the improvement from the simple Lab frame quantity A@(Lab) is moder- 
ate. The mean value of A¢(best) is 0.28, while it is 0.41 for Ag(rec, i) and 0.50 for 
A@(Lab) respectively. 

At jet level, we have to consider the possibility that the selected momenta do not 
match with parton level momenta. In Fig. B{left) and Fig. B(middle), we compare 
parton level and jet level distributions in ¢(rec)-¢(true) plain. In the left plot, we 
plot ¢(rec) with $(true) at parton level, and in the middle plot ¢(rec) is calculated for 


= J4 = 


number of events 


phi(best)—phi(true) 


Figure 7: The difference between reconstructed azimuthal angle difference and (true) at 
parton level. The solid line shows the difference from the best solution, the thick dashed 
histogram plots both of the reconstructed solutions, and the thin dotted line is the difference 
to the lab frame azimuthal angle. All histograms corresponds to 20000 events before the 
cuts. 


the two highest pr jets in the events after the detector smearing. For ¢(true) ~ 0 the 
probability to mis-reconstruct the azimuthal angle is high. This is because when two 
partons go collinear, the correlation between parton momenta and jet momenta is 
not good. In Fig H, we compare the distribution of (rec) at Point A (solid line) with 
the distribution of the events without spin correlation (a thick dashed histogram) 
and the parton level distribution without spin correlation (a thin dashed histogram). 
All distributions have a peak at ¢ = 7, but the distribution with spin correlation is 
most strongly peaked at (rec) ~ 7 as expected. 

The distribution of ¢(rec) is significantly altered by the kinematical cuts to 
reduce the background. To explain the effect of the cuts, we show the correlation 
between ¢(true) and mrp in Fig. H. In Fig. H left (middle) we show the distribution 
with (without) spin correlation. Here, mr is calculated for the two highest pr jets as 
the visible objects of the events, and Miest = Mpy. We plot both of the reconstructed 
azimuthal angle solutions in the plot. The expected endpoint of the m2 distribution 
is 600 GeV. For the two jets + missing Er signature, it is known that requiring a lower 
limit on the m2 value is an efficient cut to remove the standard model background 
BJ]. In the region ¢ ~ O(7), mrz tends to be large (small). At Point A the number of 
events at ree ~ 0 is small compared with the distribution without spin correlation. 

In Fig. P (right), we show the ¢(rec) distribution with various mrz cuts. Solid 
and dotted lines with a sharp cutoff around ¢(rec) < 2 correspond to mro(2 jet) > 
400 GeV, where mr2(2 jet) is the mr calculated with the two highest pr jets as the 
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Figure 8: Left: parton level ¢(rec) vs. (true) distribution at Point A. Middle; jet 
level d(rec) vs. (true) distribution at Point A. Right (rec) distributions at jet level 
(solid histogram), jet level without spin correlation (thick dashed histogram), parton level 
without spin correlation (thin dashed histogram). See the explanation in the text for 
details. Each plot corresponds to 2 x 104 events before the cuts. 


visible objects. About 3600 events remain after the cut for spin correlated events 
and 5000 events for spin non-correlated events. Significant suppression is observed 
for @(rec) < 2 for the events with spin correlation. The sharp edge at @ ~ 2 appears 
because the events above ¢(true) = 2 for mrz > 400 GeV are kinematically forbidden. 

The other two histograms with significant tails beyond ¢(rec) = 2 are for the 
events with mp2 >400 GeV, but the mrs is calculated for all jets in the events using 
inclusive mr2 defined in |B4, B5j. Here the jets and leptons in an event are clustered 
into two visible objects to obtain mrs. The number of events remaining after the cut 
is ~ 5300 for spin correlated events and ~ 6500 for spin non-correlated events. The 
tail events are the QQ events with significant initial state radiation. The events with 
mr2(2 jet) < 400 GeV contaminate the ¢(rec) > 2 region due to the initial state 
radiation. 

The number of events with (rec) > 2 is roughly the same between spin corre- 
lated and non-correlated samples, because ISR and the leading jets from Q decay 
are not spin correlated. This is encouraging because the events in this region may 
be regarded as control samples to determine the total number of produced events. 
Together with the number of events with ¢(rec) < 2, we may be able to determine 
how much azimuthal angle correlation exists in the events. To judge the usefulness 
of the method, it is necessary to conduct a more careful study of the dependence on 
more general production and decay patterns, and the effects of other cuts such as 
missing Er. 


5. Discussion and conclusion 


In the coming years, the LHC will accumulate more events, and the signature of 
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Figure 9: The ¢(rec) vs mr distributions with and without spin correlation (left and 
middle plots, respectively) Right: The (rec) distribution with spin correlation (solid his- 
togram) and without spin correlation (dashed histogram). Each plots generates 2 x 10* 
events before the cuts. 


supersymmetry, namely an excess of events with Ær, may be found. However, there 
are other models with similar signatures but different spin structures. In this paper 
we discuss a way to measure the spin dependence of SUSY-like signatures, especially 
those that might be observed in the early stages of the ATLAS and CMS experiments. 


We focus on the spin correlation in the 2 high pr jets + missing Er +X signature 
at the LHC. The signature would arise from a pp + QQ process when 1) the model 
has a colored quark partner Q that directly couples to quarks. 2) The colored 
particle dominantly decays into a lighter parity odd particle x and a quark q. The 
x are gauginos or Higgsinos in SUSY or heavy gauge partners in same spin partner 
(LHT or UED like) models. 3) The mass satisfies mg > my. The x may further 
decay into other particles involving the lightest “parity odd” particles. In that case, 
the jets from Q decays are prominent among the jets and initial state radiation 
[BG]. Depending on the spin of Q and its polarization, the two jets show non-trivial 
spin correlations, both in the polar angle dependence of Q — qx decay and in the 
azimuthal angle correlation of the two quarks. The effect can be seen in ETmiss VS 
Mer distributions, and the reconstructed azimuthal angle correlation using MAOS 
momentum. 

The expected spin correlation depends strongly on the production process. In 
this paper we have studied the production process based on the t channel exchange 
of the heavy gauge bosons expected in same spin partner (LHT or UED like) mod- 
els. In those models a parity odd quark Q couples directly to a quark q, and the 
interaction among Q, q and parity odd heavy gauge bosons By and Gy could be 
chiral. We have shown that the helicity structure of the amplitude is very sensitive 
to the mass and the interaction of the heavy gauge bosons. The enhanced helicity 
h = 0 component plays an especially important role. As an example, we show the 
distributions of uu => UU“), The polar angle dependence comes from overall 
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polarization, while azimuthal angle correlation comes from interference between spin 
flipping By exchange and spin conserving heavier Gy exchange. The azimuthal an- 
gle correlation is visible because the amplitude is forwardly peaked due to the By 
exchange contribution. 

Due to the azimuthal angle correlation of the events, the number of events may 
be suppressed near the mz (2 jet) end points. The high mr region is the signal 
region because the SM background is smaller, therefore spin correlations affect the 
estimation of the total production cross section. We note that there are phase space 
boundaries of the events of pp > QQ — uu + ETmiss. Some events can lie outside 
the phase space due to ISR, and they are less affected by the spin correlations. We 
note that such events may be actively used to estimate the total cross section of the 
events. 

In this paper, we do not study the corresponding distribution in SUSY models 
in detail. In supersymmetric models we do not expect spin correlations because 
the partner of the quark has spin 0. This leads to a high ETmiss/Mer, a hard mro 
end point, and a flatter reconstructed azimuthal angle correlation for two jets. Of 
course there are other important differences. the squark squark pair production is 
generally small in SUSY-like models, and gluino squark co-production is dominant. 
The gluinos may decay into squark and quark, and the ISR of the process involving 
the gluino is larger, so selecting the correct jets coming from squark decay may not 
be straightforward. However, it is shown in Ref. that the two highest pr jets in 
the events are likely from g —> xq decay. The inclusive m7. distributions with high 
pr jets show a sharp end point for the model parameter mz < mj. This means the 
two jets from squark decays are the dominant part of the event activity. Therefore 
we believe the feature discussed in this paper would be useful for SUSY studies as 
well. 

For extracting the physics behind the signatures at the LHC, using fully spin 
correlated amplitudes for production and decay is essential. To reproduce the struc- 
ture of the amplitude, the amplitude of qq > QQ > qqByBy + X should be fully 
calculated. In this paper we calculated the full amplitude using Madgraph [25], but 
other amplitude calculators or generators such as Comphep/Calchep [B] and Her- 
wig++ [B8] should be able to reproduce the effect. On the other hand, when we 
interface the two to two process to PYTHIA |R], the amplitude does not have the 
proper spin correlation, and therefore it is not appropriate for this study. 
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